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ABSTRACT 

An alternative to dark energy as an explanation for the present phase of accelerated 
expansion of the Universe is that the Friedmann equation is modified, e.g. by extra 
dimensional gravity, on large scales. We explore a natural parametrization of a gen- 
eral modified Friedmann equation, and find that the present supernova type la and 
cosmic microwave background data prefer a correction of the form 1 /H to the Fried- 
mann equation over a cosmological constant. We also explore the constraints that can 
be expected in the future, and find that there are good prospects for distinguishing 
this model from the standard cosmological constant to very high significance if one 
combines supernova data with a precise measurement of the matter density. 
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1 INTRODUCTION 

There is mounting evidence that we are living in a uni- 
verse dominated by a dark energy component, acting as 
a source of gravitational repulsion causing late-time accel- 
eration of the expansion rate. Early hints came from the 
classical test of using the magnitude-redshift relationship 
with galaxies as standard candles (Solheim 1966), but the 
reality of cosmic acceleration was not taken seriously until 
the magnitude-redshift relationship was measured recently 
using high-redshift supernovae type la (SNIa) (Riess et al. 
1998; Perlmutter et al. 1999). Cosmic acceleration requires 
a contribution to the energy density with negative pressure, 
the simplest possibility being a cosmological constant. In- 
dependent evidence for a non-standard contribution to the 
energy budget of the universe comes from e.g. the combina- 
tion of the cosmic microwave background (CMB) and the 
large-scale structure (LSS) of the Universe: the position of 
the first peak in the CMB is consistent with the universe 
having zero spatial curvature, which means that the energy 
density is equal to the critical density. Since observations 
of the LSS show that the contribution of standard sources 
of energy density, whether luminous or dark, is only a frac- 
tion of the critical density, an extra, unknown component 
is needed to account for the spatial flatness of the Universe 
(Efstathiou et al. 2002; Tegmark et al. 2003). 

The simplest explanation for the present accelerated 
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phase of expansion is to re-introduce Einstein's cosmolog- 
ical constant, A. The resulting model with A, baryons, radi- 
ation, and cold dark matter (CDM) is consistent with all 
large-scale cosmological observations like the anisotropies 
in the CMB radiation and the power spectrum of galax- 
ies (Tegmark et al. 2003). However, the value of A implied 
by the observations is tiny compared to the value inferred 
from the fact that it quantifies the energy of the vacuum in 
quantum field theory. Faced with this problem, the popular 
choice is to set A to zero and invoke a new component to 
explain the acceleration. This does not solve the problem of 
the smallness of A; if it is indeed equal to zero, one still needs 
to understand the physical mechanism behind this. Never- 
theless, one may hope that A = may be easier to explain 
than a small, but non-zero A. The question is then what the 
unknown component driving the accelerated phase of expan- 
sion is. One needs to introduce a component with negative 
pressure, and this can be done e.g. by invoking a slowly 
evolving scalar field (Wetterich 1988; Peebles & Ratra 1988; 
Ratra & Peebles 1988), or a negative-pressure fluid, e.g. a 
Chaplygin gas (Kamenshchik, Moschella & Pasquier 2001; 
Bilic, Tupper & Viollier 2001). A negative-pressure fluid, 
however, can be problematic due to its fluctuations. Fluctu- 
ations of the unknown fluid can grow very rapidly and hence 
LSS surveys can place st ri ct constraints on such models, see 
e. g. lBean fc Dord j2003l) ; ISandvik et alj (I2002T) . In order to 
circumvent this, it is sometimes assumed that the new fluid 
does not fluctuate on the scales of interest and its existence 
is visible only through modified background evolution. 

A different point of view, which we follow here, is that 
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the late time acceleration is not due to an unknown compo- 
nent but rather that the Friedmann equation is modified on 
large scales, e.g. due to extra dimensional physics. 

The structure of this paper is as follows. In section 2 
we motivate the modified Friedmann equation we discuss in 
later sections. In section 3 consider fits to SNIa data, and 
in section 4 we extend the fits to include CMB data from 
the Wilkinson Microwave Anisotropy Probe (WMAP). In 
section 5 we discuss future prospects for constraining the 
form of the modified Friedmann equation using SNIa data, 
and section 6 contains our conclusions. 



2 MODIFIED FRIEDMANN EQUATION 

We will consider a modification to the Friedmann equatio n 
with no curvature in the spirit of iDvali fc Turner! (2003), 
where they consider 



= n m (i + z) 3 + (i-a 1 



(0 



(i) 



where a is a parameter restricted to be less than 2 from 
Big Bang Nucleosynthesis (BBN) considerations. One can 
motivate Friedmann equations of this form as general 
parametrizations of the leading effect arising from modified 
gravity theories. As an example, consider a simple, single 
extra-dimensional model (Dvali, Gabadadze & Porrati 2000; 
Deffayet 2001; Deffayet, Dvali & Gabadadze 2002). The ef- 
fective, low-energy action is given by 



(2) 



S = V 1 y d A xdyy r ^)TL + J d 4 x v fg(Mp l R + Csi 

where Mpj = l/8irG, G is Newton's gravitational constant, 
g^ B is the 5-dimensional metric (A,B = 0,1,2,3,4), y is 
the extra spatial coordinate, TZ is the 5-dimensional Ricci 
scalar, g is the trace of the 4-dimensional metric, R is the 
4-dimensional Ricci scalar, and Csm is the Lagrangian of the 
fields in the Standard Model. The first term in Eq. is the 
bulk 5-dimensional Einstein action, and the second term is 
the 4-dimensional Einstein action localized on the brane at 
y = 0. The induced metric on the brane is given by g tlv (x) = 
glS{x,y — 0). The quantity r c is the new parameter of the 
theory, and is the crossover scale which sets the scale for 
the transition from 4-dimensional to 5-dimensional gravity. 
For the maximally symmetric Friedmann-Robertson- Walker 
ansatz 



ds 2 5 = f(y,H)dsl - dy 2 



(3) 



where ds\ is the 4-dimensional maximally symmetric metric, 
and H is the 4-dimensional Hubble parameter, one gets a 
modified Friedmann equation on the brane of the form 



rr2 . H 8-rrGpn 
ti zh — = 



(4) 



where p m is the matter density on the brane. This is of- 
ten called the Friedmann equation of DGP gravity (Dvali, 
Gabadadze & Porrati 2000). 

Inspired by the above example, one can also consider a 
generalized Friedmann equation 



where instead of modifying the matter content we consider 
modifications of gravity by having an arbitrary function /. 
Now assume that there is a critical scale, H c , at which 
modifications start to have an effect. Such a scale will be 
close to the present Hubble parameter. At early times, when 
H 2> H c , we know e.g. from nucleosynthesis constraints that 



f(H) r 
H c /H: 



H . In general we can then expand / in terms of 



(6) 



As long as ff > H c , evolution is standard and therefore 
co = 1 and terms with n < must vanish. Hence 



H 



1 + 



(1+*) 



(7) 



Non-standard effects start to have an effect at late times i.e 
when H ~ H c - Expanding the sum gives 



H 



H c (H c y 



= flsn»(i+z) 



(8) 



from which we see that the cosmological constant is the sec- 
ond order correction to the Friedmann equation. The first 
order correction corresponds to the DGP model. 

Generally, the nth order correction for a flat universe is 

H 



(0 



fi m (l + z) 



'(#r 



(9) 



f(H)=Hfa m (l + zf 
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which is of the same form as (0 with a = 2 — n. 

In general, one can consider constraints on the differ- 
ent coefficients c„ with n — 1,2,..., much like is done in 
parametrizing dark energy (Alam et al. 2003). In this pa- 
per, we will for simplicity only consider a single term that 
is assumed to be the leading correction. The power of the 
correction is allowed to be arbitrary i.e. we do not restrict it 
to discrete values. This approach gives information on what 
the leading order term is, and gives an idea how well one 
constrain the different terms in the expansion with current 
and future data. 



3 CONSTRAINTS FROM SUPERNOVAE 
TYPE IA 

The first test any model attempting to explain the accel- 
erated universe must pass is, of course, the SNIa data. We 
will in the fo l lowing use the sample of 194 SNIa presented in 
iBarris et alJ i2004h . The parameters we fit to the data are 
ri m and a. We consider values < fi m < 1 and — 30 < a < 2. 
The upper limit on a comes from the limits on the amount of 
energy density present at the epoch of BBN (Dvali & Turner 
2003). The Hubble parameter h is also involved in the fits, 
but is of little interest here, and we marginalize over it. The 
fit to the supernova data involves the luminosity distance 
g(l = c(l + z) J Q Z dz/H(z), and we obtain H(z) for given fi m 
and q by solving Eq. Jy with a Newton- Raphson algorithm. 
The minimum \ 2 f° r the model was 195.7 for 191 degrees 
of freedom, with the best-fitting parameters f2 m — 0.56, 
q = —14.8. The two-parameter confidence contours for f2 m 
and a are shown in Fig. It is clear that a is very weakly 
constrained by the supernova data, and can seemingly be- 
come arbitrarily large and negative. We next turn to the 
CMB data to see if they can provide firmer constraints. 
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Figure 1. 68, 95 and 99% confidence contours for f2 m and a, re- 
sulting from fitting the collection of supernova data and marginal- 
izing over h. 



Figure 2. 68, 95 and 99% confidence contours for f2 m and a, 
resulting from fitting the collection of supernova data, and adding 
the constraint from the CMB shift parameter. 



4 FITS TO CMB DATA 

Before embarking on a full fit to the CMB data, we first 
consider the much simpler approach of just fitting the so- 
called CMB shift parameter 



TZ = \/Q rn H r(zd e c), 



(10) 



where r(z) = J* dz/(H(z)/Ho) is the comoving distance 
in a flat universe, and Zdec is the redshift at decoupling. 
The shift parameter describes the shift in the CMB angu- 
lar power spectrum when the cosmological parameters are 
varied (Bond, Efstathiou & Tegmark 1997; Melchiorri et al. 
2002; Odman et al. 2003). From WMAP, 2 dcc = 1088±2, and 
ftobs = 1.716 ±0.062 (Spergel et al. 2003). Adding this con- 
straint to the supernova fit is now straightforward (Wang 
& Mukherjee 2004). For each model, we compute TZ from 
Eq. (O and add xk = (ft- ftob s ) 2 /cr^, where a n = 0.062 
to the x 2 f° r the supernova data. The resulting confidence 
contours in the ilm-a plane are shown in fig. El From the fig- 
ure it is evident that adding the constraint coming from the 
shift parameter leads to a somewhat less degenerate range 
of a, especially when considering the la contour. Still, even 
within the la limits, a can be as small as —20. Such a small 
value corresponds to a universe that becomes dominated 
by the non-standard terms in the Friedmann equation very 
quickly after H < Hq. 

The shift parameter is a useful tool for constraining cos- 
mological models quickly and very easily using information 
from the CMB. However, it does not constrain cosmological 
parameters too well. For example, in the model considered 
in this paper, considering only the shift parameter gives a 
large degeneracy along the a axis, much like in the SNIa fit. 
Hence, in order to further constrain the parameter plane, one 
needs to consider the full shape of the CMB power spectrum. 

When fitting the CMB data with a model based on an 
extra dimensional model, one should in principle start from 



the full extra dimensional theory. This can be problematic 
due to the bulk-brane interactions and hence we will in the 
following use a simplified approach where we solve the stan- 
dard 4-d perturbed Boltzmann and fluid equations, but with 
the background evolution given by Eq. ©. For convenience, 
we will parametrize the effect of the extra term in the Fried- 
mann equation by a dark energy component with an effective 
equation of state w(a). As long as the extra component does 
not fluctuate, i.e. it only has an effect on the background 
evolution, such a parametrization is equivalent to modifying 
the Friedmann equation. 

The effective equation of state is derived as follows. A 
dark energy component with equation of state w(a) has a 
density which varies with the scale factor a according to 

p D (a) = poDexp |-3^" ^-[1 + w(a')]j , (11) 

so for a flat universe, Qo = 1 — fi m , the standard Friedmann 
equation is 

(^) a = ^+(l-0„)™p{-3/"f:il + t (a')]}.(12) 

Comparing this to eq. Q, we see that we must have 

ex P {-3 J" ^l^a')]} = d3) 

for the two expressions to match. By taking the natural log- 
arithm and differentiating with respect to a on both sides, 
we get 



_l + io(o) d\n(H/H ) 

3 — = Q- 



a da 
and since da = adlna, we can write this as 

, , 1 ad\n{H/H ) 

w(a) = —1 5 — ■ -. 

v ' 3 dlna 



(14) 



(15) 
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Since a = (1 + z) , we can also write 



—— = -(l + z ) — , 
a ma dz 



giving 



w(z) 



-1 + 



a (1 + 2) d(H/H ) 



dz 



(16) 



(17) 



3 ff/Ho 

Going back to eq. Q , we can differentiate with respect to z 
and get 



H_ ^ d{H/H ) 



dz 



+ 3fi m (l + z) 



d(H/H ) 
dz 



so that 



r/z 



3fi m (l + z) 



-a(l-On)(-£- 



(18) 



(19) 



This saves us the trouble of taking numerical derivatives in 
practical work, since w(z) now can by expressed in terms of 
H(z)/H as 

aOn(l + zf 



w(z) = -1 + 



(20) 



So, instead of modifying the Friedmann equation directly, we 
can consider a standard Friedmann equation, Eq. I)12|l . with 
a fluid whose equation of state is given by (I20H . Note that 
this is an exact reformulation of the background evolution. 
We have checked that fitting the supernova data with this 
reformulation results in confidence contours in agreement 
with those in Fig. 

For fitting the CMB TT power spectrum (Hinshaw et 
al. 2003; Kogut et al. 2003) we use the likelihood code pro- 
vided by the WMAP team 1 (Verde et al. 2003). The CMB 
power spectra is computed by using CMBFAST code (Sel- 
jak & Zaldarriaga 1996), version 4.5. 1 2 . For each point in 
the parameter space (Qm,a,h), we calculate the CMB TT 
power spectrum keeping the amplitude of the fluctuations a 
free parameter by finding the best-fitting amplitude for each 
set of parameters. In other words, we only fit the shape and 
not the amplitude of the power spectrum. In calculating the 
CMB power spectrum we use Qb = 0.044 (so that we vary 
the density of cold dark matter, Q c ) and ignore reioniza- 
tion effects. Parameters Q m G [0, 1] and a € [1.5, —10] have 
uniform priors and are chosen to cover the most interest- 
ing range of parameters. For the Hubble parameter h, we 
have explored different priors: a uniform prior h £ [0.5, 1.0] 
and a Gaussian prior based on the HST Key Project value 
h = 0.72 + 0.08 (Freedman et al. 2001). We marginalize over 
h in making all the plots. The choice of prior makes little dif- 
ference after marginalization since most of the weight comes 
from h ~ 0.72. Here we show results with a Gaussian prior 
but confidence contours with a uniform prior are essentially 
identical. 

The WMAP TT power spectrum constraints are shown 
in fig. |3] From the figure we see that, as expected, having 
more information from the power spectrum than just the 
shift parameter helps tighten the constraints significantly. 




Figure 3. CMB constraints (68, 95, and 99 %) from the WMAP 
TT power spectrum 



The minimum value of a within the la contour is now only 
about —3. Comparing to fig. where only SNIa are used, 
it is obvious that CMB provides a much tighter constraint 
to the model in question than just the SNIa results. This 
demonstrates how if one restricts oneself to a single cosmo- 
logical probe in studying models with non-standard back- 
ground evolution, the CMB can be a good first choice instead 
of supernovae. If the non-standard model involves new cos- 
mological fluctuating fluids, then a simple check is provided 
by considering LSS observations which can be effective in 
constraining fluids with a non-zero sound speed. 

The combined fit CMB+SNIa is shown in fig.0] Adding 
the SNIa data further relieves the degeneracy along the a 
axis due to the fact even though both fig. Q and are both 
somewhat degenerate along a, the region of degeneracy cor- 
respond to different values of fl m . 



5 FUTURE SUPERNOVA DATA 

We have found that the present supernova data cannot put 
significant constraint on the parameter a, and that at the 
present time the most stringent constraints comes from the 
combination of CMB data with supernovae. An interesting 
question is how well one can do with future SNIa surveys. 
The planned Dark Energy Probe 3 / Supernova Acceleration 
Probe 4 is expected to observe about 2000 supernovae type 
la per year out to a redshift of z ~ 1.7 (Aldering et al. 
2002), and this should improve the power of this probe to 
constrain dark energy models considerably. We will in the 
following simulate data sets of this type, following the ap- 
proach of Saini, Weller & Bridle (2004), and consider how 



1 http://lambda.gsfc.nasa.gov/ 

2 http://www.cmbfast.org/ 



3 http:/ /universe. gsfc.nasa.gov/program/darkenergy.html 

4 http://snap.lbl.gov/ 
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Figure 4. Constraints (68, 95, and 99 %) from combining the 
WMAP TT power spectrum with supernova Type la data. 
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Figure 5. Constraints on f2 m and a from simulated data based 
on Q m = 0.3, a = 1. 



well they can constrain the the model under investigation in 
this paper. 

Empirically, SNIa are very good standard candles with 
a small dispersion in apparent magnitude <T mag = 0.15, and 
there is no indication of redshift evolution. The apparent 
magnitude is related to the luminosity distance through 



m(z) = M + 51og£> L (z), 



(21) 



where M = M + 5\og[(c/H ) Mpc" 1 ] + 25. The quantity 
Mo is the absolute magnitude of type la supernovae, and 
Dl(z) = Hodh(c)/c is the Hubble constant free luminosity 
distance. The combination of the absolute magnitude and 
the Hubble constant, M, can be calibrated by low redshift 
supernovae (Hamuy et al. 1993; Perlmutter et al. 1999). The 
dispersion in the magnitude, cr mag , is related to the deviation 
in the distance, a, by 

a - lnl0 ^ OK 

d L (z) ~ 5 mag ' [ ' 

In our simulated data sets, we assume that the errors in the 
luminosity distance are Gaussian and given by Eq. (1221 . We 
neglect systematic errors. Furthermore, we assume that the 
supernovae type la are uniformly distributed and bin them 
in 50 redshift bins, giving a relative error in the luminos- 
ity distance in each bin of ~ 1%. We do not add noise to 
the simulated cJl, and hence our results give the ensemble 
average of the parameters we fit to the simulated data sets. 

First, we simulate a data set based on a model with 
Q m ~ 0.3, a = 1. In Fig.|S]we show the confidence contours 
resulting from fitting Q m and a to this simulated data set. 
The constraints which can be derived from a data set of this 
quality are seen to be considerably tighter than those derived 
from the presently available data. However, there is still a 
notable degeneracy between f2 m and a, so that without a 
tight prior on ri m , one cannot distinguish between, e.g., a = 
1 and a = 0.5. A similar situation occurs for simulated data 
based on fi m = 0.3, a = —1, shown in Fig. |S| For this 



a -1.5 




Figure 6. Constraints on Q m and a from simulated data based 
on Qm = 0.3, a = —1. 



case, we also show in Fig. |7| the marginalized, normalized 
probability distributions for a for three different priors on 
S7 m : a uniform prior < fi m < 0.5, a Gaussian prior Q m — 
0.30 ± 0.02, and a Gaussian prior Q m = 0.300 ± 0.005. Note 
that only with the last choice of prior, where S7 m is assumed 
to be known to within 1.7 %, does one get a really tight 
constraint on a, but that even in the case of a uniform prior 
one can exclude a cosmological constant (which corresponds 
to a = 0) at 99.9 % confidence. 
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history (Lue, Scoccimarro & Starkman 2004; Multamaki, 
Gaztafiaga & Manera 2003). A more detailed analysis of 
linear and non-linear growth is left for future work. 

Finally, we note that the term we have considered is 
only the leading order correction to the standard Friedmann 
equation. The present data indicate that the first order cor- 
rection is (Ho / H) but other terms can also play a role. It is 
an interesting question whether a combination of data sets 
of the quality we can expect in the future can constrain the 
number of correction terms and their form. 
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Figure 7. Marginalized probability distribution for a with three 
different choices of priors on £7 m : uniform < Qm < 0.5 (solid 
line), Gaussian f2 m = 0.30 ± 0.02 (dashed line), and Gaussian 
n m = 0.300 ± 0.005 (dot-dashed line). 



6 CONCLUSIONS 

In this work we have studied observational constraints on 
a modified Friedmann equation that mimics dark energy in 
the universe. Modifications of the type we have considered 
here may occur naturally in models with large extra dimen- 
sions, and provide an attractive alternative to introducing an 
unknown component with negative pressure. We have found 
that the combination of the magnitude-redshift relationship 
derived from supernovae type la and the WMAP TT power 
spectrum constrain the exponent of the extra term in the 
modified Friedmann equation, (H/Ho) a , to be around — 1. 
From the point of the expanded Friedmann equation, this 
corresponds to n = 3, which is the next term after the cos- 
mological constant in the expansion (|7J . In particular, a cos- 
mological constant, corresponding to a — or n = 2 seems 
to be disfavoured for the studied parameter space. Further- 
more, the first order correction n = 1 corresponding to the 
DGP model appears to be strongly disfavored over the n = 3 
term. 

We have also used simulated data sets of the type one 
can expect from future satellite-based supernova surveys to 
estimate the accuracy with which one can hope to constrain 
the corrections to the Friedmann equation. Tight constraints 
on a can be expected if the matter density parameter f2 m is 
known accurately from other observations. But even without 
any priors on Q m , one can rule out a cosmological constant at 
high significance if the true universe is described by a = — 1. 

An obvious omission in this work is that we have not 
considered constraints coming from LSS. Since there is no 
extra negative pressure fluid in the model, one does not ex- 
pect large deviations in the matter power spectrum. Fur- 
thermore, the background evolution follows the standard 
behaviour until very recently, which suggests that linear 
growth will be standard for the most of the expansion 
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